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Abstract 

The phase space S x Z for a particle on a circle is considered. Displacement operators in this phase 
space are introduced and their properties are studied. Wigner and Weyl functions in this context 
are also considered and their physical interpretation and properties are discussed. All results are 
compared and contrasted with the corresponding ones for the harmonic oscillator in the Rx R phase 
space. 

1 Introduction 

Since the work of Wigner ^ and Moyal 0, phase space methods have been used extensively in quantum 
mechanics. A lot of this work is for the harmonic oscillator where both the position and momentum take 
values in the real line R and the phase space is the plane R x R. There has also been work on finite 
quantum systems |oJ 0] , where both the position and momentum take values in Zn (the integers module 
N) and the phase space is the lattice Zjv x Zn. The purpose of this paper is to study phase space methods 
for quantum particles on a circle. In this case the position takes values on a circle S and the momentum 
take discrete values in Z (the integers times a factor). In this case the phase space is S x Z. We note 
that in any area where there is Fourier transform involved (eg in signal processing) , the phase space can 
be R x R or Zn x Zn or S x Z in the sense that where one of the variables takes values in R or Zn or 
S 1 , the 'dual variable' take values in R or Zn or Z, correspondingly. 

Quantum mechanics on a circle is the simplest example of quantum mechanics in a non-trivial topology 
and has been studied extensively in the literature [SJE]. Physical applications include Aharonov-Bohm 
phenomena [7], mesoscopic Aharonov- Bohm rings, Floquet-Bloch wavefunctions in solid state systems, 
etc. 

In section 2 we introduce the basic formulism for position and momentum states and operators, taking 
into account the non-trivial topology of our system (described by the winding number). In section 3 we 



introduce displacement operators and study their properties. In section 4 we define Wigner and Weyl 
functions. We show that the properties of the displacement operators lead to analogous properties for 
the Wigner and Weyl functions. In section 5 we discuss an example based on a Theta wavefunction 
(which is the analogue in a circle, of a Gaussian wavefunction in a real line). Numerical examples for the 
corresponding Wigner and Weyl functions are discussed. We conclude in section 6 with a discussion of 
our results. 



2 Position and momentum states 

An electric charge is moving on a circle parameterized by the variable x. The winding number w x of x 
is defined as the integer part of the x/(2nr) (for negative x it is the integer part of the x/(2nr) minus 
1). Let r be the radius of the circle. A magnetostatic flux 4> is threading the circle in the perpendicular 
direction. The wavefunction R(x) obeys the quasi-periodic boundary condition (in units Kb = Ti = c = 1) 

R(x + 2nr) = R(x) exp(ie^), (1) 

Similar functions also appear in solid state Physics (Bloch functions). They are normalizable within each 
period 

-±-J 2 J r \R{ x) ?dx = l, (2) 
Since R(x) is a quasi-periodic function, it can be written as the following Fourier expansion: 



R(x) = R N exp(ip N x); p N = — + - ; a = (3) 

* J T ATT 



2n 

N=—od 



Rn = ~ — / exp(— ipNx)R(x)dx. (4) 



The inverse Fourier transform gives: 

2irr 

R(x) and Rn can be respectively considered as the position and momentum representations of the state 
\R). So Eq© and Eq(@J can be written as: 

oo 

(x\R) = ^2 (pn\R) exp(ip N x); (5) 

(6) 



r2nr 



1 

(Pn\R) = J (x\R) exp(-ip N x)dx. 
Let \x), \pn) be position and momentum eigenstates, correspondingly. Then 



\x) = y~] cxp(-ip N x)\p N }] \p N ) = - — / dxexp(ip N x)\x) (7) 

N=-oo J0 

(x\y) = (2irr)6 [x - y + 2Ttr(w y - w x )\ exp[-i2ira(w y - w x )); (pm\pn) = 5 M N- (8) 
It is easily seen that 

\x + 2m , w) — exp(— i2i:aw)\x) 1 (9) 
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The completeness can can be expressed as: 

i-2nr 



N=-oo 

Position and momentum operators are denned as 



— / \x)(x\dx= \pn)(pn\ = 1- (10) 

2lTr Jo 



x=- — / x\x)(x\dx, p = V pn\pn)(pn\- (11) 

We note that a different definition of x that involves integration from r to r + 2nr leads to: 

x T = x + Tl T ; II T = { \x)(x\>lx (12) 
Jo 

The x T is different from x by the projection operator II r . In the special case that r = 27rrw where w is 
an integer, n2 7rrtu = 2-irrwl. It is easily seen that: 

x\x) =(x- 2irrw x )\x), p\p N ) = Pn\Pn)- (13) 

3 Displacements and parity 

We define displacement operators as 

D(a, if) = exp ex P exp (~* a P) ' (^) 



D(a, if)|a;) = exp 



zif / . cr 
r 



(* + ?) 
) 



+ (15) 



icxK \ 

D{a,K)\p N ) = exp ( — — ) exp(-iap N )\pN+i<)- (16) 



It is easily seen that 



2r 



/if/3 Ma\ 

D(a, K)D(0, M) =D(a + f3,K + M) exp ( - —J (17) 



D(a + 27rr«7, if) = (-1)^ cxp(-i27rcru;)L>(a, if), (18) 

DHa,K) = D(-a,-K), (19) 

where w is an integer (the winding number). For later purposes we note that the D(a, if ) exp (^r-) is 
periodic in a. The period is 2irr if if is even and 4irr if if is odd number: 

i(a + 2nr)a 



D(a + 2irr, if) exp 
We also define the parity operator as: 

r2irr 



(-l) K e X p(^)D(a,K), (20) 



U = ^- \a){-a\ew( — )da= J2 \p- N )(p N \. (21) 
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The 'flux factor' exp (^r-) has been included in the definition so that the integrand \ct)(— a\ exp f— ) 
is periodic. The parity operator obeys the relations: 

U = Ul Ul = J; (22) 

The flux breaks the parity symmetry. The parity operator acting on the state |f>— jv) (which has momen- 
tum ~ N+<7 ) gives the state \pn) (which has momentum N ~^ a ); i.e., it is a parity with momentum origin 
a jr. The parity operator acting on the state \x) gives the state | — x) exp (— 2j f 2 -)- The 'flux factor' 
exp (— ^t 2 ) 'corrects' the quasi-periodicity caused by the flux into periodicity. 

We note here that in the harmonic oscillator case the displacement operators D(z) = exp(za^ — z*a) 
obey the important relations |S] 



dz R D(z) 


= V2n\p = 


>= 


-—I 


(23) 


d Zl D(z) 


= V2tt\x = 


>" 


ZR | 


(24) 


d 2 z , s 
2^ 


= u . 






(25) 



where z = Zu + izi and £/o is here the harmonic oscillator parity operator (defined in |S]). These relations 

are intimately related with the marginal properties of the Wigner and Weyl functions. 

Motivated by this we study here similar relations in our context of quantum mechanics on a circle. 

We first define the function 

1 f 27r sin(7rx) 
A(x) = — / exp(i{3x)d(3 = exp(— iirx) = exp(— «7ra;)sinc(a;). (26) 

27T J 7TX 

This is the sinc-function (used extensively in areas like digital signal processing) with a phase factor. For 
integers M, N 

A(M -N) = S(M, N) (27) 
where 8 is Kronecker's delta. Using Eas l|15|l . (|16(l can prove that: 

oo 

£ D(a,K) = \^)(~\ ) (28) 

K=-oo 

J- r r D(a,K)e X , (™) da = J ^>^' K = ^ (29 ) 

27:rJ o \ r J (En=- x \pn + k)(pn\A(-§-N), K = 2M + 1. 

We have explained earlier why the integrations involve the 'flux factor' exp (^r-)- It is natural to ask 

what is the result in Ea l29fl if we do not include the 'flux factor'. We can prove 



— / D(a,K)da= £ \p N+K ){p N \A (-- - N - a) . 

•' 1\T— _~, V / 



i-2-Kr 

* ^ / ;\ \ 

(30) 

N = ~oo 

In Ea 1)29(1 . integration from 2irr to 47rr, gives the same result for even K and the same result but with 
opposite sign for odd K: 

1 r 4 "" , r „ I |p-m)(Pm|, K = 2M; 



D(a,K)exp[ da = < (31) 

27rr7 2- -E^-oolPM>W|A(-f -JV), K = 2M+1. 



This is related to the fact (EqJ2U) that for odd K the D(a, K) exp (^) is 'anti-periodic' with 
'period' 2nr (it is periodic with period 47rr). Therefore if we integrate from to 4-7rr we get the same 
result for even K and zero for odd K. 

Combining Ecis (|21|l . (|29JI we easily show that: 

]T T r D(a 1 K)exp(—)da = Uo (32) 

^ 2nr J \ r J 

K—even 

It is natural to ask what is the result in Ea i|32[l if we sum over both even and odd integers. We can prove 

Din /V|,>Y)> ( — ) An ■ r,.4. V" V" \n, I \ ( — \f ■ V -i- ) 



Yu 2^ / D(a,K)exp(—\da = U a + £ ^ \pn+2M-i){pn\A I -M - N + - ) ; (33) 



7\ — — 00 

/>27rr 



K=-oo u v 7 M=-oo N=-oo 

The displaced parity operator is 

U(a, K) = D(a, K)U D*(a, K) = exp f £>(2a, 2iT)[/ = cxp (-^"J U D(-2a, -2K). (34) 
U(a,K) is a periodic function of a with the period irr: 

U(a + irr,K) =U(a,K). (35) 

We next show that: 

U(a,K) = -(\a)(a\+\a + Trr)(a + Trr\) , (36) 

— J U(a,K)da = \ PK )(p K \ (37) 

V— / C/(a,*Qda=l (38) 
' 27rr /n 

Eg (I36fl is proved if we use Ea (|28l) . 134|l in conjunction with the fact that 

D(a,^) = -(|-)(--| + |-+7rr){--+7rr|J . (39) 

if=even 

For completeness we also mention that 

D(a,if) = -(|-){--|-|-+7rr)(--+7rr|J . (40) 

if=odd 

Ea H4U|l can be proved using Ea (|15|l . We note that in Ea l|36|l we have projectors to the diametrically 
opposite position states |a) and \a + nr). This is related to the fact that the displaced parity operator 
on the left hand side of Ea (|36[) . has period nr. Ea l|37|l can be proved using EqJ55J). Summation over K 
in EqlO gives Eq (|35|l . 

The displacement operators arc related to the displaced parity operators through a two-dimensional 
Fourier transform. Multiplying the left and rig ht hand sides of EqlgU with D(a, K) and [D(a, K)]^ 
correspondingly, we can prove: 

' ■ >- '- r .Kf3-Ma 



M—even 



Af3=U{a,K) (41) 
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4 Wigner and Weyl functions 

The Weyl and Wigner functions can be defined in terms of the displacement and parity operator corre- 
spondingly, as: 

W(a,K) = Tr[pD(a,K)], (42) 
W(x,p N ) = Tr[pU(x,N)], (43) 

Using the fact that the density matrix is Hcrmitian, we easily show that the Wigner function is real. The 
Weyl function is in general complex. We can easily show the following formulas (which can also be used 
as alternative definitions): 

1 f 27vr a a ( K \ 

W(a,K) = — (x--\p\x+- exp [i— x ) dx, (44 
2nr ,/ 2 2 \ r J 



^2 (Pn-k\p\Pn) exp 



AT=-oo 

rl-rtr 



(45) 



W(x,pn) — - — / (x — a\p\x + a) exp(z2apAr)da, (46) 
2vrr J 

= X! (pn+k\p\pn-k) exp ( — j . (47) 

K=-oo ^ 7 ' 

The Wigner function describes the pseudo-probability of finding the particle in phase space, in a way 
consistent with quantum mechanics and the uncertainty principle. The Weyl function is equal to the 
overlap of the displaced state with the original state. In this sense, the a, K are position and momentum 
increments. The Weyl function can be understood as a generalised correlation function. If we have the 
wavefunction R(x) in order to find the correlation we displace it into R(x + a) and take the integral 
of R(x)R(x + a). In the Weyl function we perform a more general displacement in phase space i.e., a 
displacement in both position and momentum. Therefore the correlation is a special case of the Weyl 
function with a = (or K = 0). We note that the momentum takes values [N + o~)/r and depends on 
a. In contrast the momentum increments appearing on the Weyl function take values K/r and do not 
depend on a. 

The Wigner function is related to the Weyl function through a two-dimensional Fourier transform 

W(x,pn) = / > W(a, K) exp ( —i — x ) exp (iap n) da. (48) 

27rr Ir, \ r I 

This can be proved using Ea l29|) or Ea l49l) . This result is intimately connected to the fact that the 
displacement operator and displaced parity operator are related to each other through a two-dimensional 
Fourier transform fEa lf^ ). 

Using Eas l|18|) . l|35|) we easily see that the Wigner function is periodic function of x with period 7rr; 
and the Weyl function is quasi-periodic with 'period' 2irr: 

W(a + 2nrw,K) = (-l) Kw ex.p{-i2%aw)W {a, K); W(x + Trr,p N ) = W(x,p N ). (49) 



(> 



The Wigner and Weyl functions depend on the 'magnetic flux' a although for simplicity in the notation 
we have not shown this dependence explicitly. They obey the relations 

W a {a,K)=W <J+ i{a,K); W„ + i(x,p N ) = W„(x,p N+1 ). (50) 

This is related to the fact that the momentum is pn = and as we go from a to a + 1 the momentum 
Pn is relabelled as pjv+i. 

The properties of the displacement and parity operator that we proved above can be translated into 
properties for the Wigner and Weyl functions. Starting with the Weyl function we use Eas l|28|l . (|29l) to 
prove that: 

oo 

£ W(a,K) = (-|Hf); (51) 

K=-oo 

if'^il^M,. JW*-* K = ™> (52) 

2 " J « V ' ) I E».- K (Pivl*» + riA(-f -AT), K = 211+1. 



We can make here a similar comment to that made for Ea (|29(l . In Ea (|52|) . integration from 2irr to Anr, 
gives the same result for even K and the same result but with opposite sign for odd K . Therefore if we 
integrate from to Airr we get the same result for even K and zero for odd K . We can also show that 

W{Q,Q) = l, W(a,K) = W*(-a,-K), \W(a,K)\<l. (53) 

For the Wigner function we use Eas (|36|) . H37[l . (|38|) to prove: 

— J W{x,p N )dx = (pn\p\pn), (54) 

^ W(x,p N ) = -((x\p\x) +(x + Trr\p\x + Ttr)), (55) 



2 

JV=-oo 



'0 N 

We can also prove that 



, „2ixr oo 

/ J2 W(x, PN )dx = 1. (56) 
Z7Tr Jo „_ „ 



Tr(pip 2 ) = t, — / Yl W pi( x 'PN)W P2 (x,p N )dx, (57) 

JU N=-oo 

5 Example 

As an example, we consider a 'free particle' descibed with the Hamiltonian 

H = f. (58) 

We note that the particle feels a vector potential through the quasi-periodic boundary condition. 

We assume that at t = the wavefunction is a 'Gaussian on a circle' i.e., a Theta function. In order 
to explain this we introduce the Zak transform [|] on a function S(y) on a real line, defined as 

oo 

R(x, a) =JV S(y = x+ 2irrw) exp(-i27rcrw), (59) 



W — — OQ 
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where N is used to normalize the function R(x,a) according to EqJ2J). If S(y;A) is a Gaussian wave- 
function 

S(y; A) = 7T" 1 / 4 exp f-?L + ^2A ■ y - AAj^j (60) 
where A — Ar + iAj, then the Zak transform is 



R(x,a; A) = Af ■ 7r~ 1/4 exp f-y + V2Ax - AAr) 9 3 -ttct + mr{x - V2A); i2nr 2 (61) 
where Q^[u;t] is Theta function ^Uji defined as 

oo 

©3[u;r]= exp(i7TTn 2 + ilnu). (62) 



n— — oo 



We have taken r = 1, A = 1 and calculated the Wigner and Weyl functions. Numerical results are shown 
in Figs 1-3. 

In Fig.l we present a contour plot of the Wigner function W(x, N + a) as a function of x and for all 
values of a. For a given a the Wigner function is defined only on the discrete values of the momenta 
Pn = (N + a)/r and the parallel black lines in the figure show the case a — 0.1. As expected from Ea H49|) . 
the Wigner function is periodic function of x with period it. Also it has been explained in Ea (|50|) that 
the plotted values of W(x,pi) for a = 1 represent the W(x,p2) for a = 0. 

In Fig. 2 we present the Wigner function as a function of a; at a particular momentum p$ — a with 
a = 0.1. This is really an appropriate slice of the three-dimensional version of Fig.l. 

In Fig. 3 we present a contour plot of the absolute value of the Weyl function as a function of a for 
K = 1 and for all values of a. The black line in the figure shows the case a = 0.1. 

We calculate the time evolution of this system. This is easily done in the momentum representation 

as 

( PN \R(t)) = cM*tH)(p N \R{0)) (63) 

We need the state of the system at t = in the momentum representation. If S(p) is the Fourier transform 
of S(x): 

r+oc 

S{p) = / S(x) exp(-ipx)dx 1 (64) 



then using EqQJ, Ea (|59|l and Ea i|tj4|) we can prove the relation: 



Rn = ^S{ PN ). (65) 
Zirr 



Since the Fourier transform of Gaussian wavefunction Eg 1601) is: 

S( P ; A) = V27T 1 / 4 exp (- ?L - iV2A ■ p + AAj j . ( W, i 

we obtain: 

(p N \R) = -jL— gxpf-ZlL-iy/zA-pN + AAi ) . (<>7) 



s 



Inserting Ea (|57|) into we get 



A/V 1 / 4 
(p N \R(t)) = — = — exp 



it ) p N 2 - iV2A ■ p N + AAi 



(68) 



From that we have calculated for t — 1 the corresponding Wigner and Weyl functions to the previous 
ones. Results are shown in Figs 4-6. The results show the time evolution of the system in the language 
of Wigner and weyl functions. 



6 Discussion 

Quantum mechanics on a circle has attracted a lot of attention in the literature. It uses topological ideas 
in the context of quantum mechanics. In this paper we have studied phase space methods in this context. 
We have introduced displacement operators and studied their properties. There is clearly some analogy 
with the harmonic oscillator displacement operator properties (given in Eas (|23|) . H24[l . (|25[) % ). but there are 
also considerable differences. In particular the flux a plays an important role in the properties of the 
displacement operators for particles on a circle. Our main results here are Eas (l28ll . (|29[) . (|32|) . 

We also introduced Wigner and Weyl funcitons and discussed their physical interpreration and their 
properties, which are direct consequence of the displacement operator properties. A numerical example 
for the Theta wavefunction of Ea (|61|l (the analogue on a circle of the Gaussian function on the real line) 
has been discussed. 

The results can be used in the context of Aharonov-Bohm devices, Floquet-Bloch solid state systems, 
coherent states on a circle quantum maps and classical and quantum chaos |12| etc. Fractional 
Fourier transforms in this context have been discussed in |13) . 
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Figure 3: A contour plot of the absolute value of the Weyl function at t = for K = 1, as a function of 
a and a. The black line in the figure shows the case a = 0.1 
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Figure 6: A contour plot of the absolute value of the Weyl function at t = 1 for K = 1, as a function of 
a and a. The black line in the figure shows the case a — 0.1 
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